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Abstract 

Graphene's honeycomb lattice structure is quite remarkable in the sense that it leads, in the long 
wavelength limit, to a massless Dirac equation description of nonrelativistic quasiparticles associ- 
ated with electrons and holes present in the two dimensional crystallite. In the case of cold bosonic 
atoms trapped in a honeycomb optical lattice, Haddad and Carr (2009) have recently shown, by 
taking into account binary contact interactions, that the dynamics of these Bose-Einstein conden- 
sates is governed by a nonlinear Dirac equation (NLDE). In this paper, we study exact stationary 
solutions of such a NLDE. After proving that the energy eigenvalues are real, we show that the sum 
of orbital angular momentum and pseudospin angular momentum normal to the crystal commutes 
with the nonlinear Hamiltonian whenever magnitudes of the pseudospin components do not depend 
on the polar angle <j). We obtain some exact stationary and localized solutions of the NLDE. 



I. INTRODUCTION 



Graphene, peeled off from bulk graphite by means of micromechanical cleavage technique, 
is a zero-gap semiconductor [1,2]. Being a monolayer of carbon atoms arranged in a planar 
hexagonal lattice structure in which charge carriers travel distances that are thousand times 
the lattice spacing without being subjected to any scattering, graphene helped in estab- 
lishing that stable two dimensional crystals do exist, contrary to what one expects from 
thermodynamic instability arguments put forward by Landau, Peierls and Mermin in the 
past [3,4]. 

For wavelengths much larger than the lattice size, charged quasiparticles in the honey- 
comb lattice satisfy a massless Dirac equation [5-7], 



dtp _, 
z/t— = —ihvF(a.V)tp (1) 




where, 



dx ^ dy 
a = ia x + ja y 

and vf ~ 0.01 c represents the Fermi speed. The components ip A and if)B are bi-component 
spinors describing spin-half charge carriers, while subscripts A and B refer to the two in- 
equivalent sites A and B of the sublattice forming the hexagonal structure, representing 
pseudospin degrees of freedom. The associated Clifford algebra, 

yy + 7"r" = 2 rT, /*, f = 0, 1, 2 

is satisfied by, 

7° = <7z, 7 1 = <7z(*x, 7 2 = <*z<* V ■ 

Relativistic quantum features like zitterbewegung and Klein paradox, which so far have been 
confined to the domain of quantum electrodynamics, make their appearance in the context 
of subrelativistic charge carriers in graphene because of eq.(l) [8]. 

In an interesting recent paper, Haddard and Carr have studied weakly interacting bosonic 
atoms anchored to a planar honeycomb optical lattice [9] . The quantum many body theory 



of this system is described by a Hamiltonian, 

H = J J&(^^V 2 + V(r?)$dxdy + \J J V$t$$ dx dy ( 2 ) 



where the field of identical atoms each having mass m is represented by a complex scalar field 
operator $(r, t). The potential V(f) encapsulates interactions of atoms with electromagnetic 
fields from laser beams responsible for pinning the atoms to a two dimensional optical lattice. 
The second term in eq.(2) arises because of binary contact interaction between atoms, A 
being a coupling constant proportional to the s-wave scattering length. 

As a result of both the planar honeycomb lattice structure as well as the presence of a 
contact interaction term in eq.(2), the dynamics of Bose-Einstein condensate (BEC) in this 
case is described by a nonlinear Dirac equation, 

ih^t = if(V# = -ihc s (a.V)ip + M($)ij) (3) 

where c s ~ few cm s _1 is an effective BEC sound speed in the planar lattice, and ip represents 
the state of BEC, 



if, 



/ tp A (x,y,t) 



V 



M(j/j) = U 



Mx,y,t) 

with ipA and if)B now being scalar wavefunctions representing spatio-temporal and pseudospin 
degrees of freedom for the bosonic atoms, while M(ip) embodies the nonlinearity arising out 
of s-wave scatterings, 

v \M 2 ; 
U being the interaction energy [9]. 

From eq.(3), one can straight away obtain the conservation of probability current, 
where, 

P = = V^A + 

and, 

J = c s ifj^aip 

Because of eq.(4), one may normalize the BEC wavefunction, 

J Jpdxdy = J J(\ip A \ 2 + \ij B \ 2 )dxdy = l (5) 
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II. DISPERSION RELATION, INNER PRODUCT AND ENERGY 



In order to study dispersion relation in NLDE, we substitute the following plane wave 
form in eq.(3), 

(J^ e i(k x x+k y y~u>t) 
(J Be i(k x x+k y y-wt) 



1> 



that entails, 



and, 



(u - ^r\C B \ 2 )C B = c s (k x + ik y )C A (6a) 
U 



\lo - ^\C A \ 2 )C A = c s (k x - ik y )C B (66) 
Assuming box-normalization for ip on a square of side length L, one obtains from eq.(5), 

\C A \ 2 + \C B \ 2 = ^ (7) 

Making use of eqs. (6a), (6b) and (7), we arrive at, 



where, 



Reality of \C A \ 2 requires, 



A* = u '-t>M + k;). (9) 



hu ^ 1 h 2 L 2 A 2 

T*u? + -ip- (10) 



Therefore, for U > 0, eq.(10) implies that if, 

U h 2 L 2 A 2 , , 

plane wave solutions of eq.(3), corresponding to 'massive modes' with tiny mass \J^-0~, are 
possible. 

If the inner product of two BEC states ip and \ is defined in the usual manner as, 

{^,x) = / / ^Xdxdy (12) 

one can trivially verify that, 

(M(^)x,^) = (x,M(^) (13) 
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which entails, 

(H(4>) X ,i>) = (x,H(4>W) (14) 

since —ih(a.W) is self-adjoint. 

Therefore, if -0^ is an energy eigenstate so that, 

H(iJ>e)iI>e = EipE (15) 

then substituting ip E in place of \ and ip in eq.(14), leads to the result that E is real. 
However, if 

H(tpi)^i = Eii>i , t = l,2 

then, in general, ipi and ip2 are not mutually orthogonal for non-zero U. 

(It must be pointed out here that technically we should not be using terms like eigenvalues 
and eigenstates since we are dealing with a nonlinear operator H{ip). But since these terms 
have entered the common parlance of quantum theory, we will continue to use them.) 

Suppose U = so that the Hamiltonian appearing in eq.(3) is just — iTic s {a.V). Then, 
it is interesting to observe that the orbital angular momentum normal to the planar lattice 
does not commute with the Hamiltonian, 

[L z , -iha.V] = h 2 [o x — - a y —J (16) 

Now, since, 

/ d d \ 

[a z ,-iha.V] = -2h[a x - VyTr) (17) 



dy dx 



we can get using eqs.(16) and (17), 



[L z + ^a z ,-iM.V]=0. (18) 

suggesting that L z + |<r 2 plays the role of total angular momentum. But this is puzzling as a 
appearing in eq.(3) reflects only pseudospin degrees of freedom. In the context of graphene, 
using much more extensive arguments, Mecklenburg and Regan have recently made a case 
for identifying pseudospin with intrinsic spin angular momentum [10]. Although eq.(3) 
describes a bosonic system, the pseudospin degrees arising out of sites A and B of a two 
dimensional honeycomb lattice appear to induce a half-integral spin angular momentum for 
the BEC when U = 0. 
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III. LOCALIZED STATIONARY SOLUTIONS 



In this section, we look for energy eigenfunctions of the NLDE whose magnitudes fall with 
radial distance. It is convenient to switch to polar coordinates (r, (p), so that a stationary 
state solution tp E (r : (p : t) of the NLDE has the form, 

Vfc = e-^fc (p) 

where, 

ip(r, (p) -- 

satisfies the time independent NLDE, 

E^^ Cs e-^-^) + U\ M ^ (19) 



Mr, </>) 



EipB = -^.e*(^T + % -^f) + U\1>b\ 2 1>b , (20) 



' dtp A i dtp A ' 
dr r d(p 

E being the energy. 

In polar coordinates, the orbital angular momentum is simply L z = —ih-^. Therefore, 

[L z ,M(^)]=0 

only when both \ip A \ 2 and \iPb\ 2 are independent of the angle (p. Since [a z , M(ip)\ vanishes 
for all ip, in order that L z + |cr z commutes with H(ip), the two components of ip must take 
the following form, 

i> A (r,<j>) = R A (r)e iQ ^ (21) 

and, 

ifj B {r,(P) = R B {r)e ieB(r ^ (22) 

We seek solutions that are simultaneous eigenfunctions of the nonlinear Hamiltonian and 
the total angular momentum normal to the optical lattice. So, making use of eqs.(21) and 
(22) in eqs.(19) and (20) we find, 

& A (r, ( p) = (l- 1 -) ( p+fA(r) (23) 
e B (r,0) = (Z + ^ + / B (r) (24) 
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for / = 0, 1, 2..., while the radial parts satisfy, 

ER A = -ihc s e^-M ( ^± + lRB d JjL + ^L (l + M + UF p (25 ) 

V dr ar r 2 / 

£i2 fl = -ihc.e^°-M(^± + iRa^ - ^(l - \)) + (26) 

When Ra = Rb, using the reality of energy E, one can deduce from eqs. (25) and (26) 
the following class of exact solutions, 

R A = R B = Kr- 1 ' 2 (27) 

f B (r) = f A (r) + Q (28) 

where, 

<?.*,-' (-^), (29) 

with being the normalization constant. 
If I = 0, we get, 

f A (r) = ^-r 1 - [Er - UK 2 In r}+K x (30) 

tlCg 

and, 

/B(r) = /A(r) + n7r (31) 

where K\ is an arbitrary real constant. 

For other integral values of I, eqs. (25)-(27) demand vanishing of energy eigenvalue along 
with, 

UK 2 

f A (r) = —j—cosQhir + K 2 , (31) 
where Q is given by eq.(29) and K 2 is a real constant. 

If we normalize the wavefunction over a circular region of radius L, the constant K gets 
fixed to, 

K 1 



y/AnL 

To summarize our results, we have obtained exact solutions of a nonlinear Dirac equation 
given by eq.(3) for the case \i/ja\ = \4>b\ that correspond to definite energy and angular 
momentum normal to the plane of the optical lattice. The position probability density for 
these solutions decrease as r _1 . While the s-wave solution can have any value of energy, 
solutions with 1 — 1,2.. are static in nature since E — 0. 



What is interesting is that the planar honeycomb lattice not only leads to a massless 
nonlinear Dirac equation description of the BEC but also makes the pseudospin degrees of 
freedom appear as half- integral spin degrees in a purely bosonic system (also see [10]). 
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